Highly damped quasinormal modes of Kerr black holes 



m 
o 
o 

(N 

o 



E. Berti^^^, V. Cardoso^, K.D. Kokkotas^ and H. Onozawa"^ 
^ Department of Physics, Aristotle University of Thessaloniki, Thessaloniki 54124, Greece 
^ Centra Multidisciplinar de Astrofisica - CENTRA, Departamento de Fisica, 
Institute Superior Tecnico, Av. Rovisco Pais 1, 1049-001 Lishoa, Portugal 
^ San Diego Wireless Center, Texas Instruments, 
5505 Morehouse Dr., San Diego, CA, 92121, USA 
* McDonnell Center for the Space Sciences, Department of Physics, 
Washington University, St. Louis, Missouri 63130, USA * 
(Dated: February 1, 2008) 

Motivated by recent suggestions that highly damped black hole quasinormal modes (QNM's) may 
provide a link between classical general relativity and quantum gravity, we present an extensive 
computation of highly damped QNM's of Kerr black holes. We perform the computation using 
two independent numerical codes based on Leaver's continued fraction method. We do not limit 
our attention to gravitational modes, thus filling some gaps in the existing literature. As already 
observed in [isll . the frequency of gravitational modes with I = m = 2 tends to lor — 2Q, Q being 
the angular velocity of the black hole horizon. We show that, if Hod's conjecture is valid, this 
asymptotic behaviour is related to reversible black hole transformations. Other highly damped 
modes with m > that we computed do not show a similar behaviour. The real part of modes 

Owith I — 2 and m < seems to asymptotically approach a constant value ljh ~ ~mvj, zu ~ 0.12 
being (almost) independent of a. For any perturbing field, trajectories in the complex plane of 
' QNM's with m — show a spiralling behaviour, similar to the one observed for Reissner-Nordstrom 

(RN) black holes. Finally, for any perturbing field, the asymptotic separation in the imaginary part 
^SJ . of consecutive modes with m > is given by 2hTh (Th being the black hole temperature). We 

^ ' conjecture that for all values of I and m > there is an infinity of modes tending to the critical 

CO ' frequency for superradiance {ujr = m) in the extremal limit. Finally, we study in some detail 

modes branching off the so-called "algebraically special frequency" of Schwarzschild black holes. 
■ For the first time we find numerically that QNM multiplets emerge from the algebraically special 

Schwarzschild modes, confirming a recent speculation. 

o ■ 
cn 

^ ■ 1. INTRODUCTION 

The study of linearized perturbations of black hole solutions in general relativity has a long history |lj. The 
Oh! development of the relevant formalism, initially motivated by the need of a formal proof of black hole stability, gave 
birth to a whole new research field. A major role in this field has been played by the concept of quasinormal modes 
(QNM's): oscillations having purely ingoing wave conditions at the black hole horizon and purely outgoing wave 
conditions at infinity. These modes determine the late-time evolution of perturbing fields in the black hole exterior. 
Numerical simulations of stellar collapse and black hole collisions in the "full" (non-linearized) theory have shown that 
in the final stage of such processes ( "ringdown" ) QNM's dominate the black hole response to any kind of perturbation. 
Since their frequencies are uniquely determined by the black hole parameters (mass, charge and angular momentum), 
QNM's are likely to play a major role in the nascent field of gravitational wave astronomy, providing unique means 
to "identify" black holes 0. 

An early attempt at relating QNM's to the Hawking radiation was carried out by York jsj. More recently Hod made 
an interesting proposal to infer quantum properties of black holes from their classical oscillation spectrum . It was 
suggested many years ago by Bekenstein Q that in a quantum theory of gravity the surface area of a black hole (which 
by the Bekenstein-Hawking formula is nothing but its entropy) should have a discrete spectrum. The eigenvalues of 
this spectrum are likely to be uniformly spaced. Hod observed that the real parts of the asymptotic (highly damped) 
quasinormal frequencies of a Schwarzschild black hole of mass M, as numerically computed by Nollert 0] and later 
by Andersson 'T'l , can be written as 

ujb = ThI'D-3, (1) 

where we have used units such that c = G = 1 and Th is the black hole's Hawking temperature. He then exploited 
Bohr's correspondence principle, requiring that "transition frequencies at large quantum numbers should equal classical 
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oscillation frequencies" , to infer that variations in the black hole mass induced by quantum processes should be given 
by 

AM = huR. (2) 

Finally, he used the first law of black hole thermodynamics to deduce the spacing in the area spectrum for a 
Schwarzschild black hole. Remarkably, in this quantum gravity context relevant modes are those which damp in- 
finitely fast, do not significantly contribute to the gravitational wave signal, and are therefore typically ignored in 
studies of gravitational radiation. Following Hod's suggestion, Dreyer recently used a similar argument to fix a free 
parameter (the so-called Barbero-Immirzi parameter) appearing in Loop Quantum Gravity j^] . Supposing that transi- 
tions of a quantum black hole are characterized by the appearance or disappearance of a puncture with lowest possible 
spin jmim Dreyer found that Loop Quantum Gravity gives a correct prediction for the Bekenstein-Hawking entropy 
if jmin = 1, consequently fixing the Barbero-Immirzi parameter. 

When Hod made his original proposal, formula was merely a curious numerical coincidence. Kunstatter Q 
suggested that a similar relation may hold also for multidimensional black holes. Since these early speculations, a 
full formalism for non-rotating black hole perturbations in higher dimensions has been developed [lOj, and different 
calculations have now shown that formula Q holds exactly for scalar and gravitational perturbations of nonrotating 
black holes in any dimension [T2i. .13. ,14. 15. l&j. Furthermore, Birmingham et al. have recently given intriguing 
hints corroborating the correspondence suggested by Hod focusing attention on 2-1-1 dimensional Bahados- 
Teitelboim-Zanelli (BTZ) black holes In this case the QNM frequencies (which belong to two "families" ) can be 
obtained analytically, and their real parts are independent of the mode damping. They showed that the identification 
of the fundamental quanta of black hole mass and angular momentum with the real part of the QNM frequencies 
leads to the correct quantum behaviour of the asymptotic symmetry algebra, and thus of the dual conformal field 
theory. 

In light of these exciting new results, Hod's conjecture seems to be a very promising candidate to shed light on 
quantum properties of black holes. However, it is natural to ask whether the conjecture applies to more general 
(charged and/or rotating) black holes. If asymptotic frequencies for "generic" black holes depend (as they do) on the 
hole's charge, angular momentum, or on the presence of a cosmological constant, should Hod's proposal be modified in 
some way? And how does the correct modification look like? The hint for an answer necessarily comes from analytical 
or numerical calculations of highly damped QNM's for charged and rotating black holes, or for black holes in non- 
asymptotically flat spacetimcs. Some calculations in this direction have now been performed, revealing unexpected 
and puzzling features JJ,, 19, 20, 21, 22, 23, 24, 25, 26,, 27J. 

In particular, the technique originally developed by Nollert to study highly damped modes of Schwarzschild black 
holes has recently been extended to the RN case [l^, showing that highly damped RN QNM's show a peculiar spiralling 
behaviour in the complex-w plane as the black hole charge is increased. Independently, Motl and Neitzke obtained 
an analytic formula for the asymptotic frequencies of scalar and electromagnetic-gravitational perturbations of a RN 
black hole whose predictions show an excellent agreement (at least for large values of the charge) with the numerical 
results [T^ . For computational convenience they fixed their units in a somewhat unconventional way: they introduced 
a parameter k related to the black hole charge and mass by Q/M = 2y/k/{l + k), so that f3 = 47r/ (1 — fc) = 1/Th is 
the inverse black hole Hawking temperature and Pi = —k'^fi is the inverse Hawking temperature of the inner horizon. 
Their result is an implicit formula for the asymptotic QNM frequencies, 

e'^" 2 + 'de-f^"^ = 0, (3) 

which has recently been confirmed by independent calculations |24|. However, its interpretation in terms of the 
suggested correspondence is still unclear. Asymptotic quasinormal frequencies of a charged black hole, according to 
formula depend not only on the black hole's Hawking temperature, but also on the Hawking temperature of 
the (causally disconnected) inner horizon. Perhaps more worrying is the fact that the asymptotic formula does not 
yield the correct Schwarzschild limit as the black hole charge Q tends to zero. The mathematical reason for this 
behaviour has been discussed in 0,0] • A calculation of higher-order corrections in uij may explain the observed 
disagreement: indeed, as we shall see, the numerical study of Kerr modes we present in this paper seems to support 
this expectation. Finally and most importantly, it is not at all clear which are the implications of the generally non- 
periodic behaviour of asymptotic RN modes for the Hod conjecture. Maybe the complicated behaviour we observe is 
an effect of the electromagnetic-gravitational coupling, and we should only consider pure gravitational perturbations 
for a first understanding of black hole quantization based on Hod's conjecture. The latter suggestion may possibly 
be ruled out on the basis of two simple observations: first of all, in the large damping limit "electromagnetic" and 
"gravitational" perturbations seem to be isospectral to each other, and isospectral to scalar perturbations as well jl2l |: 
secondly, Kerr modes with m = show a very similar spiralling behaviour, which is clearly not due to any form of 
electromagnetic-gravitational coupling. 
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The available numerical calculations for highly damped modes of black holes in non-asymptotically flat spacetimes 
are as puzzling as those for RN black holes in flat spacetime. Cardoso and Lemos [23| have studied the asymptotic 
spectrum of Schwarzschild black holes in a de Sitter background. They found that, when the black hole radius is 
comparable to the cosmological radius, the asymptotic spectrum depends not only on the hole's parameters, but 
also on the angular separation index I. Formula does not depend on dimensionality and gives the same limit for 
"scalar" and "gravitational" modes (loosely using the standard four dimensional terminology; see 0, 0| for a more 
precise formulation in higher dimensions). This "universality" seems to be lost when the cosmological constant is 
non-zero. The study carried out in ,^0] has recently been generalized to higher dimensional Schwarzschild-de Sitter 
black holes [2^ and to take into account higher-order corrections to the predicted behaviour . However the issue 
is not settled yet, and the asymptotics may be different from what was predicted in pol |. Indeed, recent numerical 
and analytical calculations [23, Hj] seem to suggest that the result presented in is only correct when the overtone 
index n satisfies nfc <C 1, where k is the surface gravity at the Schwarzschild-de Sitter black hole horizon. For higher 
overtones, the behaviour seems to be different. The problem is not completely solved yet. Numerically, it seems 
difficult to compute QNM frequencies for nk > 1 Furthermore, at present, numerical and analytical results show 
only a qualitative (but not quantitative) agreement |27| . 

Calculations of QNM's for Schwarzschild-anti-de Sitter black holes were performed in various papers (22j, showing 
that the nature of the QNM spectrum in this case is remarkably different (basically due to the "potential barrier" 
arising because of the cosmological constant, and to the changing QNM boundary conditions at infinity). Those 
calculations were recently extended to encompass asymptotic modes |23| . The basic result is that consecutive highly 
damped modes (whose real part goes to infinity as the imaginary part increases) have a uniform spacing in both the 
real and the imaginary part; this spacing is apparently independent of the kind of perturbation considered and of the 
angular separation index I. 

The aim of this paper is to study in depth the behaviour of highly damped Kerr QNM's, complementing and 
clarifying results that were presented in previous works 0, |2^ . The plan of the paper is as follows. In section [21 
we briefly introduce our numerical method. In sectional we discuss some results presented in and show a more 
comprehensive calculation of gravitational QNM's, considering generic values of m and higher multipoles (namely, 
I = S). In section 21 we display some results for scalar and electromagnetic perturbations. If our numerics for 
non-gravitational modes are indicative of the true asymptotic behaviour, the asymptotic formula which is valid for 
I = m = 2 gravitational perturbations may be very special. In section \^ we briefly summarize our results and we 
discuss the asymptotic behaviour of the modes' imaginary part. Finally, in section we turn our attention to a 
different open problem concerning Kerr perturbations. Motivated by some recent, surprising developments arising 
from the study of the branch cut in the Schwarzschild problem |30| and by older conjectures derived from analytical 
calculations of the properties of algebraically special modes [sj], we turn our attention to Kerr QNM's in the vicinity 
of the Schwarzschild algebraically special frequencies. As the black hole is set into rotation, we flnd for the flrst time 
that a QNM multiplet appears close to the algebraically special Schwarzschild modes. A summary, conclusions and 
an outlook on possible future research directions follow. 



2. NUMERICAL METHOD 

A first numerical study of Kerr QNM's was carried out many years ago by Detweiler [s^- Finding highly damped 
modes through a straightforward integration of the perturbation equations is particularly difficult even for non- 
rotating black holes '2*1 . In the Kerr case the situation is even worse, because, due to the non-spherical symmetry 
of the background, the perturbation problem does not reduce to a single ordinary differential equation for the radial 
part of the perturbations, but rather to a system of differential equations (one equation for the angular part of the 
perturbations, and a second equation for the radial part). 

A method to find the eigenfrequencies without resortin g to integrations of this system was developed by Leaver, 
and has been extensively discussed in the literature |l9l l29l l33| . In this paper we will apply exactly the same method. 
Following Leaver, we will choose units such that 2M = 1. Then the perturbation equations depend on a parameter 
s denoting the spin of the perturbing field (s = 0, —1, —2 for scalar, electromagnetic and gravitational perturbations 
respectively), on the Kerr rotation parameter a (0 < a < 1/2), and on an angular separation constant Aim- In 
the Schwarzschild limit the angular separation constant can be determined analytically, and is given by the relation 
Aim + s{s + l). 

The basic idea in Leaver's method is the following. Boundary conditions for the radial and angular equations 
translate into convergence conditions for the series expansions of the corresponding eigenfunctions. In turn, these 
convergence conditions can be expressed as two equations involving continued fractions. Finding QNM frequencies 
is a two-step procedure: for assigned values of a, £, m and w, first find the angular separation constant Aimi^^) 
looking for zeros of the angular continued fraction equation; then replace the corresponding eigenvalue into the radial 
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continued fraction equation, and look for its zeros as a function of lu. Leaver's method is relatively well convergent 
and numerically stable for highly damped modes, when compared to other techniques [s^l- We mention that an 
alternative, approximate method for finding Kerr quasinormal frequencies has recently been presented j^^j which has 
the advantage of highlighting some physical features of the problem. 

In the next sections we will use Leaver's technique to complement numerical studies of Kerr quasinormal overtones 
carried out by some of us in the past 0, . The method we use for our analysis is the one described in those 
papers. Exploring the high-damping regime necessarily requires pushing our numerics to their limits. Therefore we 
have systematically cross-checked the reliability of our results using two independent codes. As we shall see, our 
study will uncover a plethora of interesting new features. 



3. GRAVITATIONAL PERTURBATIONS 



3.1. I — m — 2 modes: a more extensive discussion 



Let us consider rotating black holes, having angular momentum per unit mass a = J/M. The black hole's (event 
and inner) horizons are given in terms of the black hole parameters by r± = M ± V — a^. The hole's temperature 
Th ~ {r+ — rJ)IA where A = 87rAfr+ is the hole's surface area, related to its entropy S by the relation S = A/4. 
Introducing the angular velocity of the horizon 17 = 47ra/A, applying the first law of black hole thermodynamics, 

AM = ThAS* + r^A J, (4) 

and assuming that the formula for the area spectrum derived for a Schwarzschild black hole still holds in this case, 
Hod conjectured that the real parts of the asymptotic frequencies for rotating black holes are given by: 

i^R — i^R = Th In 3 + mfi, (5) 

where m is the azimuthal eigenvalue of the field 0. Hod later used a systematic exploration of moderately damped 
Kerr black hole QNM's carried out a few years ago by one of us 29] to lend support to formula 0, at least for modes 
with I = m ISp . His conclusions were shown to be in contrast with the observed behaviour of modes having stronger 
damping in |l9j|: the deviations between the numerics and formula (jSJ were indeed shown to grow as the mode order 
grows (see figure 7 in p^). Hod even used equation |0J, without including the term due to variations of the black hole 
charge AQ, to conjecture that © holds for Kerr-Newman black holes as well ,4]. This second step now definitely 
appears to be a bold extrapolation. Not only does formula (O disagree with the observed numerical behaviour for 
perturbations of Kerr black holes having I = m = 2 [l^ (not to mention other values of m, as we shall see below); 
by now, analytic and numerical calculations have shown that RN QNM's have a much more rich and complicated 
behaviour 

In summary, there is now compelling evidence that the conjectured formula jSJ) must be wrong. However it turns 
out Il9j, quite surprisingly, that an extremely good fit to the numerical data for ^ = m = 2 is provided by an even 
simpler relation, not involving the black hole temperature: 

ujii = toO. (6) 

At first sight, the good fitting properties of this formula may be regarded as a coincidence. After all, this formula 
does not yield the correct Schwarzschild limit. Why should we trust it when it is only based on numerical evidence? 
A convincing argument in favour of formula (jSJ is given in figure ^ There we show the real part of modes having 
Z = m = 2 as a function of n for some selected values of a (namely, a = 0.05,0.10, ..,0.45). The convergence towards 
the limiting value lvr = 2^1 (horizontal lines in the plot) is evident. Furthermore, the convergence is much faster 
for holes spinning closer to the extremal limit, and becomes slower for black holes which are slowly rotating. The 
behaviour we observe presents interesting analogies with the asymptotic formula lj2Jl. The Schwarzschild limit r nay 
not be recovered straightforwardly as a — > 0. Some order-of- limits issues may be at work, as recently claimed in |l3| 
to justify the incorrect behaviour of formula Q as the black hole charge Q ^ 0. 

Is formula ® merely an approximation to the "true" asymptotic behaviour, for example a lowest-order expansion 
in powers of f2? To answer this questions we can try and replace © by some alternative relation. Since in the 
Schwarzschild limit equation ((BJ doesn't give the desired "In 3" behaviour, we would like a higher-order correction 
which does reproduce the non-rotating limit, while giving a good fit to the numerical data. Therefore, in addition to 
equations Q and we considered the following fitting relations: 

ujj^^ = 47rT^ln3 + mf2 = rffln3(l - fi^) +mfi, (7) 
CUR = TH\n3{l-m^n'^) + mn. (8) 
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Formula lO enforces the correct asymptotic limit at a = 0, and can be considered as an r2^-correction to Hod's 
conjectured formula ((SJ). Since numerical results suggest a dependence on mQ we also used the slight modification 
given by formula © , hoping for a better fit to our numerical data. The relative errors of the various fitting formulas 
with respect to the numerical computation for the n — 40 QNM are given in figure |21 Equation JHJl is clearly the one 
which performs better. All relations are seen to fail quite badly for small rotation rate, but this apparent failure is 
only due to the onset of the asymptotic behaviour occurring later (that is, when n > 40) for small values of a. 

We believe that the excellent fitting properties and the convergence plot, when combined together, are very good 
evidence in favour of equation ©. Maybe the impressive visual agreement between the numerics and the conjectured 
asymptotic formula ©, displayed in the left panel of figure 13 is even more convincing. Therefore, let us assume as 
a working hypothesis that equation © is the correct asymptotic formula (at least for Z = m = 2, and maybe for 
large enough a), and let us consider the consequences of such an assumption in computing the area spectrum for Kerr 
black holes. Modes having I = m may indeed be the relevant ones to make a connection with quantum gravity, as 
recently claimed in (36(|. Furthermore, the proportionality of these modes to the black hole's angular velocity fl seems 
to suggest that something "deep" is at work in this particular case. 

In the following, we will essentially repeat the calculation carried out by Abdalla et al. |^ for near-extremal 
(a M) Kerr black holes. We will argue that the conclusion of their calculation is in fact wrong, since those authors 
did not take into account the functional behaviour of ^^(a) (which was unknown when they wrote the paper), but 
rather assumed that lor — mjlM is constant in the vicinity of the extremal limit. In following the steps traced out 
in [221 we will restore for clarity all factors of M. This means, for example, that the asymptotic frequency for m > 
in the extremal limit is lor — to/2 A/. Let us also define x — a/M. The black hole inner and outer horizons are 
r± = M [1 ± (1 — x^)^/^] . The black hole temperature is 



and we recall that the black hole surface area A = SnM^ [l + (1 — a;^)^/'^] is related to its entropy S by the relation 
S = A/ A. The hole's rotational frequency is 

^^4™^^^J X ^ 

A 2Mr+ 2Ml + VT^ 

Let us now apply the first law of black hole thermodynamics and the area-entropy relation to find 

AA = ^ (AM - nAJ) . (11) 

The authors of focused on the extremal limit. They used AJ = km and AM = fujJii{x = 1) = hm/2M to deduce 
that 



A^ = Ahm 



1/2M -VL 
T 



= hmA. (12) 



where A is the area quantum. Now, the square parenthesis is undefined, since VI — > 1/2 Af when x — > 1 . Taking the 
limit X ^ 1 and keeping AM = hm/2M constant leads to 



A^87:\^l + ^^-^j ^87:, (13) 

which is the final result in |22j |. The fundamental assumption in this argument is that the asymptotic frequency is 
ujji = m/2M^ which is strictly true only for x = 1. However, one has to consider how the QNM frequency changes 
with X. What is the effect of assuming u)r — mil on the area spectrum? The calculation is exactly the same, but the 
equation AM = hm/2M is replaced by AM — hmil, and we conclude that 

AA = 0. (14) 

The area variation is zero at any black hole rotation rate a < M. At first sight, this result may look surprising. 
It is not, and it follows from fundamental properties of black holes. Indeed, we are looking at reversible black hole 
transformations. It is well known that the gain in energy AE and the gain in angular momentum A J resulting from 
a particle with negative energy —E and angular momentum —L^ arriving at the event horizon of a Kerr black hole is 
subject to the inequality 

AM > QAJ; (15) 
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see, for example, equation (352) on page 373 in [j| and the related discussion. This inequality is equivalent to the 
statement that the irreducible mass M^^r = (Mr+/2)^/^ of the black hole can only increase |33- In other words, 
by no continuous infinitesimal process involving a single Kerr black hole can the surface area of the black hole be 
decreased (Hawking's area theorem). Assuming the validity of Hod's conjecture and using the result © for 
asymptotic QNM's, we are saturating the inequality IjlSfl: we are considering a reversible process, in which the area 
(or, equivalently, the irreducible mass) is conserved. Classically, this result makes sense. Perturbations of Kerr black 
holes dying out on a vanishingly small timescale are likely to be a process for which the horizon area is an adiabatic 
invariant. Some physical processes exhibiting this feature were considered in detail in [3^ . 

What does the result H14(l mean from the point of view of area quantization? It could mean that using modes 
having Z = m in Hod's conjecture is wrong, or that we cannot use Bohr's correspondence principle to deduce the 
area spectrum for Kerr. A speculative suggestion may be to modify Bohr's correspondence principle as introduced by 
Hod. Suppose for example that we do not interpret the asymptotic frequencies as a change in mass (AM = fiLOR), 
but rather impose TAS — hujR. This is of course equivalent to Hod's original proposal when a — 0. The asymptotic 
formula would then imply, using the first law of black hole thermodynamics, that the minimum possible variation in 
mass is AM = 2mMl. 

We notice that the above arguments do not apply to strictly extremal Kerr black holes, for which a — M. In the 
extremal case the horizon area is not an adiabatic invariant |39j , and its quantization probably requires some special 
treatment. 

3.2. Modes with I = 2, I ^ m 

As discussed in the previous paragraph, we feel quite confident that the real part of modes with I ^ m = 2 
approaches the limit loh = wf2 as the mode damping tends to infinity. What about modes having I ^ ml In [l9l | 
it was shown that modes with m = show a drastically different behaviour. As the damping increases, modes show 
more and more loops. Pushing the calculation to very large imaginary parts is not easy, but the trend strongly 
suggests a spiralling asymptotic behaviour, reminiscent of RN modes. In this section we present results for the cases 
not considered in [13, concentrating on the real parts of modes with I = 2 and m — 1, — 1, — 2. 

Modes for which I = 2, m = 1 are displayed in the right panel of figure |3 They do not seem to approach the limit 
one could naively expect, that is, ujr — Q. Instead, the real part of the frequency shows a minimum as a function of 
a, and approaches the limit lur — m as a 1/2. To our knowledge, the fact that the real part of modes with I — 2 
and m = 1 approaches ujr = m= lasa^l/2 has not been observed before. In the following we will see that this 
behaviour is characteristic of QNM's due to perturbation fields having arbitrary spin, as long as m > 0. 

The real parts of modes with I — 2, m < as functions of a (for some selected values of n) are displayed in figure^l 
From the left panel, displaying the real part of modes with m = — 1, we infer an interesting conclusion: the frequencies 
tend to approach a constant (presumably a-independent) limiting value, with a convergence rate which is faster, as 
in the I = m = 2 case, for large a. The limiting value is approximately given by 0.12. A similar result holds for modes 
with I — 2, m — —2 (right panel). Once again the frequencies asymptotically approach a (roughly) constant value, 
with a convergence rate which is faster for large a. The limiting value is now approximately given by ujr — 0.24, 
about twice the value we got for m = — 1. In summary, the real part of modes with m < seems to asymptotically 
approach the limit 

lor ~ — mi37, (16) 

where vo ~ 0.12 is (to a good approximation) independent of a, at least in the extremal limit a 1/2. 

We will see below that this surprising result is quite general. It is supported by calculations of gravitational QNM's 
for different values of I, and it also holds for electromagnetic and scalar perturbations, as long as m < 0. An analytical 
derivation of this result is definitely needed. It may offer some insight on the physical interpretation of the result, and 
help explain the surprising qualitative difference in the asymptotic behaviour of modes having different values of to. 

3.3. Modes with / = 3 

Results for a few highly-damped QNM's with Z = 3, to = were shown in Those modes exhibit the usual 

"spiralling" behaviour in the complex plane as the imaginary part increases. In this paragraph we present a more 
complete calculation of modes with I = 3. Some care is required in considering the results of this section as represen- 
tative of the asymptotic behaviour. In fact the pure imaginary Schwarzschild algebraically special mode (separating 



7 



the lower QNM branch from the upper branch) is located at 

^ . ii-mi + i){i + 2) 

and can be taken as (roughly) marking the onset of the asymptotic regime. The algebraically special mode quickly 
moves downwards in the complex plane as I increases, and corresponds to an overtone index n = 41 when I = 3. 
Unfortunately we did not manage to push our numerical calculations for I = 3 to values of n larger than about 50. 
Therefore we cannot be completely sure that our results are indicative of the "true" QNM asymptotics. 

In any event, some prominent features emerge from the general behaviour of the real part of the modes, as displayed 
in the different panels of figure [3 First of all, contrary to our expectations, neither the branch of modes with m — 3 
nor the branch with m = 2 seem to approach the limit we would expect, luh = mn. These modes show a behaviour 
which is more closely reminiscent of modes having I — 2, m — 1: the modes' real part "bends" towards the zero- 
frequency axis, shows a minimum as a function of a, and tends to lur = m as a — > 1/2. If the qualitative behaviour 
of QNM's does not drastically change at larger overtone indices, we would be facing a puzzling situation. Indeed, 
gravitational modes with / = m = 2 would have a rather unique asymptotic behaviour, that would require more 
physical understanding to be motivated. 

Another prominent feature is that, whenever m > 0, there seems to be an infinity of modes approaching the limit 
ujfi = m as a 1/2. This behaviour confirms the general trend we observed for I = 2, m > 0. 

Finally, our calculations of modes with to < show, once again, that these modes tend to approach LUf; — —mw, 
where w ~ 0.12. We display, as an example, modes with I = 3 and m = —1 in the bottom right panel of figure[Sl 



4. SCALAR AND ELECTROMAGNETIC PERTURBATIONS 



The calculations we have performed for I — 3 hint at the possibility that modes with I = m = 2 are the only ones 
approaching the limit lor = mn. However, for reasons explained in the previous paragraph, carrying out numerical 
calculations in the asymptotic regime when I > 2 is very difhcult. 

This technical difficulty is a hindrance if we want to test the "uniqueness" of gravitational modes with I = m ^ 2 
by looking at gravitational modes having I > 2. An alternative idea to check this "uniqueness" is to look instead 
at perturbations due to fields having different spin and I < 2. In particular, here we show some results we obtained 
extending our calculation to scalar (s = 0)and electromagnetic (s = —1) modes. To our knowledge, results for Kerr 
scalar modes have only been published in [3a | . Some highly damped electromagnetic modes were previously computed 




4.1. Scalar modes 

In figure we show a few scalar modes with / = m = 0. As we could expect from existing calculations p^, "s^ the 
modes show the typical spiralling behaviour; the surprise here is that this spiralling behaviour sets in very quickly, 
and is particularly pronounced even if we look at the first overtone {n — 2). As the mode order grows, the number of 
spirals grows, and the centre of the spiral (corresponding to extremal Kerr holes) moves towards the pure imaginary 
axis (at least for n < 10). 

In figure [7| we show the trajectories of some scalar modes for 1 = 2. As can be seen in the top left panel, rotation 
removes the degeneracy of modes with different values of to. If we follow modes with to = we see the usual spiralling 
behaviour, essentially confirming results obtained in j35j using the Priifer method. However our numerical technique 
seems to be more accurate than the (approximate) Priifer method, and we are able to follow the modes up to larger 
values of the rotation parameter: compare the bottom right panel in our figure 7 to figure 6 in |35j |. and remember 
that their numerical values must be multiplied by a factor 2 (due to the different choice of units) . On the basis of our 
numerical results, it is quite likely that the asymptotic behaviour of scalar modes with Z = to = is described by a 
relation similar to ||2Jl. However, at present, no such relation has been derived analytically. 

In figure IHlwc show the real part of scalar modes with I = m = 1 and Z = to = 2 as a function of a, for increasing 
values of the overtone index n. In both cases modes do not show a tendency to approach the lor = mfl limit 
suggested by gravitational modes with I = m = 2. As we observed for modes with I = 3 and to > 0, their behaviour 
is rather similar to that of gravitational modes with 1 = 2 and to = 1. This may be considered further evidence that 
gravitational perturbations with I = m = 2 are, indeed, very special. 
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4.2. Electromagnetic modes 

The calculation of highly damped electromagnetic QNM's basically confirms the picture we obtained from the 
computation of scalar QNM's presented in the previous section. We show some selected results in figure El The top 
left panel shows that, for large damping, the real part of electromagnetic QNM's with / = 1 and m > shows a 
local minimum, approaching the limit lor = m as a — > 1/2. The top right panel shows that the real parts of modes 
with / = 1 and to = quickly start oscillating (that is, QNM's display spirals in the complex-w plane). Finally, the 
bottom plots show the behaviour of modes with Z = 1, to = —1 (left) and I = 2, m = —2 (right). Once again, if our 
calculations arc indicative of the asymptotic behaviour, modes seem to approach a roughly constant value uj ~ —mvo. 

5. THE ASYMPTOTIC BEHAVIOUR OF THE MODES' IMAGINARY PART 

The evidence for a universal behaviour emerging from the calculations we have presented is suggestive. For reasons 
we explained in the previous sections, in some instances we may not have reached the asymptotic regime when our 
numerical codes become unreliable. With this caution, we can still try and draw some conclusions. Our results suggest 
that, whatever the kind of perturbation (scalar, electromagnetic or gravitational) that we consider, asymptotic modes 
belong to one of three classes: 

1) Modes with to > 0: their real part probably approaches the limit — mQ only for gravitational modes with 
I = TO. Our calculation for / = to = 3 cannot be considered as a trustworthy counterexample to this prediction, 
since it is not really representative of the asymptotic regime. For other kinds of perturbations (and for m ^ I) lor 
apparently shows a minimum as a function of a. This may be a real feature of asymptotic modes, but it may as well 
be due to the asymptotic behaviour emerging only for larger values of n. To choose between the two alternatives we 
would either require better numerical methods or the development of analytical techniques. A "universal" feature is 
that, whatever the spin of the perturbing field, QNM frequencies approach the limiting value ujr = to as a ^ 1/2. 

2) Modes with to = 0: these modes show a spiraling behaviour in the complex plane, reminiscent of RN QNM's. 

3) Modes with to < 0: their real part seems to asymptotically approach a constant (or weakly a-dependent) limit 
LOji ~ — TOtu, where vj ~ 0.12, whatever the value of I and the spin of the perturbing field. Maybe this limit is not 
exactly independent of a, but on the basis of our numerical data we are quite confident that highly damped modes 
with TO < tend to a universal limit ojr ~ —mwext, where vjext has some value between 0.11 and 0.12, as a ^ 1/2. 

Another interesting result concerns the modes' imaginary part. In [T9| we observed that the following formula holds 
for gravitational modes with I = m = 2: 

^/=™=2 = 2fi + i2TrTHn. (18) 

Our numerical data show that, in general, all modes with to > have an asymptotic separation equal to 2ttTh- 
This result holds for all kinds of perturbations (scalar, electromagnetic or gravitational) we considered, as long as 
TO > 0. For TO = the imaginary part oscillates, and this beatiful, general result does not hold. It turns out that it 
doesn't hold as well for modes with to < 0. So far the analysis of our numerical data did not lead us to any conclusion 
on the asymptotic separation of modes with to < 0. This may hint at the fact that for m < our calculations are 
not yet indicative of the asymptotic regime. Therefore, some care is required in drawing conclusions on asymptotic 
modes from our results for to < 0. 



6. ALGEBRAICALLY SPECIAL MODES 

6.1. An introduction to the problem 

Algebraically special modes of Schwarzschild black holes have been studied for a long time, but only recently their 
understanding has reached a satisfactory level. Among the early studies rank those of Wald 40] and of Chandrasekhar 
[ilf. who gave the exact solution of the Regge- Wheeler, Zerilli and Teukolsky equations at the algebraically special 
frequency. The nature of the QNM boundary conditions at the Schwarzschild algebraically special frequency is 
extremely subtle, and has been studied in detail by Maassen van den Brink |3l|. Black hole oscillation modes belong 
to three categories: 

1) "standard" QNM's, which have outgoing wave boundary conditions at both sides (that is, they are outgoing 
at infinity and "outgoing into the horizon" , using Maassen van den Brink's "observer-centered definition" of the 
boundary conditions); 
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2) total transmission modes from the left (TTM^'s) are modes incoming from the left (the black hole horizon) and 
outgoing to the other side (spatial infinity); 

3) total transmission modes from the right (TTM^'s) are modes incoming from the right and outgoing to the other 
side. 

In our units, the Schwarzschild "algebraically special" frequency is given by formula H17|l . and has been traditionally 
associated with TTM's. However, when Chandrasekhar found the exact solution of the perturbation equations at the 
algebraically special frequency he did not check that these solutions satisfy TTM boundary conditions. In (3l)| it was 
shown that, in general, they do not. An important conclusion reached in 31] is that the Regge- Wheeler equation 
and the Zerilli equation (which are known to yield the same QNM spectrum, being related by a supersymmetry 
transformation) have to be treated on different footing at fi; , since the supersymmetry transformation leading to the 
proof of isospectrality is singular there. In particular, the Regge- Wheeler equation has no modes at all at fi;, while 
the Zerilli equation has both a QNM and a TTMl ■ 

Numerical calculations of algebraically special modes have yielded some puzzling results. Studying the Regge- 
Wheeler equation (that should have no QNM's at all according to Maassen van den Brink's analysis) and not the 
Zerilli equation. Leaver found a QNM which is very close, but not exactly located at, the algebraically special 
frequency. Namely, he found QNM's at frequencies llj such that 

= 0.000000 - 3.9980001, tl'^ = -0.000259 - 20.0156531. (19) 

Notice that the "special" QNM's Cl[ are such that 3? iJlj < \^2\, ^ ^^'3 > l^sl, and that the real part of Ct'^ is not 
zero. Maassen van den Brink [3ll | speculated that the numerical calculation may be inaccurate and the last three 
digits may not be significant, so that no conclusion can be drawn on the coincidence of Qi and fJJ, "if the latter does 
exist at all". 

An independent calculation was carried out by Andersson Using a phase-integral method, he found that the 
Regge- Wheeler equation has pure imaginary TTM^^'s which are very close to il/ for 2 < / < 6. He therefore suggested 
that the modes he found coincide with f2;, which would then be a TTM. Maassen van den Brink observed that, 
if all figures in the computed modes are significant, the coincidence of TTM's and QNM's is not confirmed by this 
calculation, since f2J and Qi are numerically (slightly) different. 

Onozawa showed that the Kerr mode with overtone index n — 9 tends to fi; as a 0, but suggested that 
modes approaching 17/ from the left and the right may cancel each other at a = 0, leaving only the special (TTM) 
mode. He also calculated this (TTM) special mode for Kerr black holes, solving the relevant condition that the 
Starobinsky constant should be zero and finding the angular separation constant by a continued fraction method; his 
results improved upon the accuracy of those previously obtained in [4Tl |. 

The analytical approach adopted in clarified many aspects of the problem for Schwarzschild black holes, but the 
situation concerning Kerr modes branching from the algebraically special Schwarzchild mode is still far from clear. In 
|3lj| Maassen van den Brink, using slow-rotation expansions of the perturbation equations, drew two basic conclusions 
on these modes. The first is that, for a > 0, the so-called Kerr special modes (that is, solutions to the condition 
that the Starobinsky constant should be zero '25,l4l|) are all TTM's (left or right, depending on the sign of s). The 
TTM/j's cannot survive as a — s- 0, since they do not exist in the Schwarzschild limit; this is related to the limit a — s- 
being a very tricky one. In particular, in this limit, the special Kerr mode becomes a TTM^ for s = —2; furthermore, 
the special mode and the TTM/? cancel each other for s = 2. Studying the limit a ^ in detail, Maassen van den 
Brink reached a second important conclusion: the Schwarzschild special frequency fli is a limit point for a multiplet 
of "standard" Kerr QNM's, which for small a are well approximated by 

33078176 3492608. , 

Lj = -4i maH ia^ + O(ma^) + O(a^) 20) 

700009 41177 ^ v ; -r v / v 7 

when / = 2, and by a more complicated formula - his equation (7.33) - when I > 2. None of the QNM's we numerically 
found seems to agree with the analytic prediction when the rotation rate a is small. 

Maassen van den Brink suggested (see note [46] in (23|) that QNM's corresponding to the algebraically special 
frequency with m > may have one of the following three behaviours in the Schwarzschild limit: they may merge 
with those having m < at a frequency flj such that < jfi/j (but Jfij] > \ for ^ > 4) and disappear, as suggested 
by Onozawa pol ]: they may terminate at some (finite) small a; or, finally, they may disappear towards lu — — ioo. 
Recently Maassen van den Brink et al. [sol ] put forward another alternative: studying the branch cut on the imaginary 
axis, they found that in the Schwarzschild case a pair of "unconventional damped modes" should exist. For / = 2 
these modes were identified by a fitting procedure to be located on the unphysical sheet lying behind the branch cut 
(hence the name "unconventional") at 



uj± = T0.027+ (0.0033- 4)1. 



(21) 
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An approximate analytical calculation confirmed the presence of these modes, yielding 

UJ+ ~ -0.03248 + (0.003436 - 4)i, (22) 

in reasonable agreement with H21|l . If their prediction is true, an additional QNM multiplet should emerge near Q,i as 
a increases. This multiplet "may well be due to uj± splitting (since spherical symmetry is broken) and moving through 
the negative imaginary axis as a is tuned" |30j . In the following paragraph we will show that a careful numerical 
search indeed reveals the emergence of such multiplets, but these do not seem to behave exactly as predicted in [30i] . 

6.2. Numerical search and QNM multiplets 

As we have summarized in the previous paragraph, the situation for Kerr modes branching from the algebraically 
special Schwarzschild mode is still unclear, and there are still many open questions. Is a multiplet of modes emerging 
from the algebraically special frequency when a > 0? Can QNM's be matched by the analytical prediction H20I) at 
small values of a? If a doublet does indeed appear, as recently suggested in jSOj, does it tend to the algebraically 
special frequency VI2 = — 4i as a ^ 0, does it tend to the values predicted by formula (|21|l . or does it go to some other 
limit? 

After carrying out an extensive numerical search with both our numerical codes, we have found some surprises. 
Our main new result is shown in the left panel of figure EH There we show the trajectories in the complex plane of 
QNM's with I — 2 and to > 0: a doublet of modes does indeed appear close to the algebraically special frequency. 
Both modes in the doublet tend to the usual limit {(I2 — ni) as a — > 1/2. We have tried to match these "twin" modes 
with the predictions of the analytical formula (|20|l . Unfortunately, none of the two branches we find seems to agree 
with H2()|l at small a. Our searches succeeded in finding a mode doublet only when to > 0. For to < the behaviour 
of the modes is, in a way, more conventional. For example, in the right panel of figure^] we see the ^ = 2, to = 
mode emerging from the standard algebraically special frequency ^2 and finally describing the "usual" spirals as a 
increases. 

In the top left panel of figurelTTlwe see that the real part of all modes having to > does indeed go to zero as a ^ 0, 
with an TO-dependent slope. However, the top right panel in the same figure shows that the imaginary part of the 
modes does not tend to —4 as a — > 0. Qualitatively this behaviour agrees rather well with that predicted by equation 
(|21(l . Extrapolating our numerical data to the limit a — > yields, however, slightly different numerical values; our 
extrapolated values for I = 2 are ui = (—4 — 0.10)i and oj = (—4 + 0.09)i. 

At present, we have no explanation for the appearance of the doublet only when to > 0. A confirmation of this 
behaviour comes from numerical searches we have carried out for I =^ 3, close to the algebraically special frequency 
ris. Once again, a QNM multiplet only appears when to > 0. In particular, we see the appearance of a doublet that 
behaves quite similarly to the modes shown in the left panel of figure ITUI Extrapolating the numerical data for the 
I = 3 doublet yields the values uj = (-20 - 0.19)i and w = (-20 + 0.24)i as a ^ 0. 

A more careful search near the algebraically special frequency O3 surprisingly revealed the existence of other QNM's. 
However, the additional modes we find may well be "spurious" modes due to numerical inaccuracies, since we are 
pushing our method to its limits of validity (very high dampings and very small imaginary parts). 



7. CONCLUSIONS 



In this paper we have numerically investigated the behaviour of highly damped QNM's for Kerr black holes, using 
two independent numerical codes to check the reliability of our results. Our findings do not agree with the simple 
behaviour conjectured by Hod for the real part of the frequency ^4, 36J as given in equation (jsj. We did not limit our 
attention to gravitational modes, thus filling some gaps in the existing literature. 

Our main results concerning highly damped modes can be summarized as follows. Scalar, electromagnetic and 
gravitational modes show a remarkable universality of behaviour in the high damping limit. The asymptotic behaviour 
crucially depends, for any kind of perturbation, on whether to > 0, to = or to < 0. As already observed in |l9j |. 
the frequency of gravitational modes with I — m — 2 tends to ujn = 2Q, Q being the angular velocity of the black 
hole horizon. We showed that, if Hod's conjecture is valid, this asymptotic behaviour is related to reversible black 
hole transformations, that is, transformations for which the black hole irreducible mass (and its surface area) does 
not change. 

Other (gravitational and non-gravitational) modes with to > do not show a similar asymptotic behaviour in the 
range of n allowed by our numerical method. In particular, in the high-damping limit, the real part of (gravitational 
and non-gravitational) modes with to > typically shows a minimum as a function of the rotation parameter a, and 
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then approaches the Umit lob. = m as the black hole becomes extremal. At present we cannot exclude the possibility 
that our calculations actually break down before we reach the asymptotic regime. Better numerical methods or 
analytical techniques are needed to give a final answer concerning the asymptotic behaviour of modes with m > 0. 

Hod recently used a continued-fraction argument modelled on that used in Tl'l and claimed that the asymptotic 
Kerr QNM frequency is given (for any m) by 

tjKorr ^ _^ YIttThU. (23) 

This result is obviously compatible with our calculations only for to > 0, so there is some reason to be cautious about 
Hod's derivation. Essential in his argument is a comparison of the order of magnitude of the recursion coefficients a„, 
(3n and 7n defined in equations (6), (7) and (8) of (4^. Looking at his formulas (14) and (15), it is apparent that the 
magnitudes of the a„ and 7„ recursion coefficients as n — > oo are of the same order. However, in the «„ terms are 
treated as negligible with respect to the and 7„ terms. Equation H23|) comes from imposing 7n = for n > N, where 
iV is some (large) integer. However, if 7„ = for all n > N, then it is not legitimate to say that the 7n-terms are much 
larger than the a„-terms. Neglecting the a„-terms is not correct: after imposing 7„ = 0, the expansion coefficients 
dn for large n are computed by comparing the q;„ and /3„ terms, not the f3n and 7n terms. Furthermore, if Hod's 
argument were correct, it would allow the calculation of the real part of the frequency for Schwarzschild gravitational 
perturbations. However, applying his argument to the Schwarzschild case, Hod could only derive the asymptotic 
behaviour of electromagnetic perturbations, for which the QNM frequency vanishes. Finally, a contradiction with our 
numerics would result if the asymptotic limit were reached for n 3> (a6)~^, where b = (1 — 4a^)^/^, as stated in [43l |. 
Our numerical results show that this is only valid for I = m = 2, otherwise we would clearly see convergence to the 
asymptotic behaviour already for n = 30 — 50 (at least for intermediate values of a and b). 

Recently Musiri and Siopsis showed that equation ((SJ holds in an intermediate regime, when |a;| is large but \uja\ < 1 
|44j . Their result is compatible with our calculations, and (unfortunately) it does not provide a final answer on the 
asymptotic behaviour. Concluding, despite these recent efforts, a more careful analytical analysis is needed before 
drawing any final conclusion on asymptotic Kerr QNM frequencies. 

An interesting new finding of this paper is that for all values of to > 0, and for any kind of perturbing field, there 
seems to be an infinity of modes tending to the critical frequency for superradiance, ujh = to, in the extremal limit. 
This finding generalizes a well-known analytical result by Detweiler for QNM's with I = to (32, ||35|]. It would be 
interesting to generalize Detweiler's proof, which only holds for ^ = to, to confirm our conjecture that for any to > 
there is an infinity of QNM's tending to ojr = to as a ^ 1/2. 

The real part of modes with I = 2 and negative to asymptotically approaches a value ujn ~ —mzu, w ~ 0.12 being 
(almost) independent of a. Maybe this limit is not exactly independent of a, but on the basis of our numerical data 
we feel confident that highly damped modes with to < do tend to a universal limit w/j ~ —mwext (where tUext has 
some value between 0.11 and 0.12) as a — > 1/2. This is an interesting prediction, and it would again be extremely 
useful to confirm it using analytic techniques. So far we have not been able to find any simple physical explanation for 
this limiting value. For example, we have tentatively explored a possible connection between w and the frequencies 
of marginally stable counterrotating photon orbits, but we could not find any obvious correlation between the two. 

Both for gravitational and for non-gravitational perturbations, the trajectories in the complex plane of modes with 
TO = show a spiralling behaviour, strongly reminiscent of the one observed for Reissner-Nordstrom (RN) black holes, 
and probably well approximated in the high damping limit by an equation similar to 

Last but not least, an important result concerning highly damped modes is that, for any perturbing field, the 
asymptotic separation in the imaginary part of consecutive modes with to > is given by 2ttTh {Th being the black 
hole temperature). An heuristic explanation for this fact was put forward for the Schwarzschild case in |l3|. The idea 
is as follows. Since QNM's determine the position of the poles of a Green's function on the black hole background, 
and the Euclidean black hole solution converges to a thermal circle at infinity having temperature Th, it may not be 
too surprising that the spacing in asymptotic QNM's coincides with the spacing 2'k\Th expected for a thermal Green's 
function. However, this simple relation concerning the mode spacing does not seem to hold when to < 0. Analytic 
derivations for the spacing in the QNM imaginary parts have been provided in 45] and 46] . These calculations use 
the fact that QNM's are poles in the scattering amplitude of the relevant wave equation. They are based on the Born 
approximation, and they only apply to static spacetimes. A generalization to stationary spacetimes, if possible, might 
provide an analytical confirmation of our numerical result. 

Finally, we studied in some detail modes branching from the so-called "algebraically special frequency" of 
Schwarzschild black holes. We found numerically for the first time that QNM multiplets emerge from the alge- 
braically special modes as the black hole rotation increases, confirming a recent speculation [33 . However, we found 
some quantitative disagreement with the analytical predictions in |3Cll3l| . The problem deserves further investigation. 

Hopefully our numerical results will serve as a guide in the analytical search for asymptotic QNM's of Kerr black 
holes. Although one can in principle apply Motl and Neitzke's 12] method in the present case, the Kerr geometry 
has some special features that complicate the analysis. The Tcukolsky equation describing the field's evolution no 



12 



longer has the Regge- Wheeler- Zerihi (Schrodinger-like) form; however, it can be reduced to that form by a suitable 
transformation of the radial coordinate. The main technical difficulty concerns the fact that the angular separation 
constant Aim is not given analytically in terms of ^, as it is in the Schwarzschild or RN geometry; even worse, it depends 
on the frequency w in a non linear way. Therefore, an analytical understanding of the problem must also encompass 
an understanding of the asymptotic properties of the separation constant. The scalar case is well studied, both 
analytically and numerically |47| . but a similar investigation for the electromagnetic and gravitational perturbations 
is still lacking. An idea we plan to exploit in the future is to use a numerical analysis of the angular equation as a 
guideline to find the asymptotic behaviour of Aim- Once the asymptotic behaviour of Aim is determined, the analysis 
of the radial equation may proceed along the lines traced in |12j. 
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FIG. 1: Each different symbol corresponds to the (numerically computed) value of ljr as a function of the mode index n, at 
different selected values of the rotation parameter a. The selected values of a are indicated on the right of the plot. Horizontal 
lines correspond to the predicted asymptotic frequencies 2fl at the given values of a. Convergence to the asymptotic value is 
clearly faster for larger a. In the range of n allowed by our numerical method (n < 50) convergence is not yet achieved for 
a < 0.1. 
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FIG. 2: Relative difFerence between various fit functions and numerical results for the mode with overtone index n = 40. From 
top to bottom in the legend, the lines correspond to the relative errors for formulas 0, (|SJ, Q and ||5J. 
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FIG. 3: Real part of the frequency for different modes with I — 2 and m > 0. In both panels we overplot (bold solid line) the 
prediction of formula (|5J. The left panel shows the excellent agreement between modes with I = m = 2 and the asymptotic 
formula. The right panel shows the different behaviour of modes with m = 1; these modes have a frequency that "bends" 
downwards as n increases, showing a local minimum as a function of a. In both cases, uir ^ m in the extremal limit a —> 1/2. 





FIG. 4: Real part of the first few modes with I — 2 and m < 0. Modes with m — —1 are shown in the left panel, modes with 
m = —2 in the right panel. As the mode order n increases, uir seems to approach a (roughly) constant value uiu — —mw, 
where vj ~ 0.12. Convergence to this limiting value is faster for large values of the rotation parameter a (compare figure 0. 




FIG. 5: Real parts of some modes with I 
behaviour is reminiscent of modes with I - 
LOR = —mzo (we only show modes with m 



! = 3 and different values of m (indicated in the plots). When m > 0, the observed 
= 2, m = 1 (see figure |^. Modes with m < approach a (roughly) constant value 
= —1), as they do for / = 2 (see figure |1J. 
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FIG. 7: Trajectories of a few scalar modes for 1 = 2. In the top left panel we show how rotation removes the degeneracy of 
modes with different m's, displaying three branches (corresponding to m = 2, 0, — 2) "coming out of the Schwarzschild limit" 
for the fundamental mode (n = 1). In the top right and bottom left panel we show the progressive "bending" of the trajectory 
of the m = branch for the first two overtones (n = 2, 3). Finally, in the bottom rig ht p anel we show the typical spiralling 
behaviour for a mode with m = and n = 9. This plot can be compared to figure 6 in |3fll | (notice that their scales have to be 
multiplied by two to switch to our units). The continued fraction method allows us to compute modes for larger values of a 
(and is presumably more accurate) than the Priifer method. 
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FIG. 8: Real parts of the scalar modes with i = m = 1 (left), I = m = 2 (right). The observed behaviour is reminiscent of 
figures 121 and 1^ 
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FIG. 10: The left panel shows the trajectories described in the complex-w plane by the doublets emerging close to the 

Schwarzschild algebraically special frequency (O2 = — 4z) when rn > and I — 2. Notice that the real part of modes with 
m > tends to ujr = m as a ^ 1/2. The right panel shows the spiralling trajectory of the mode with m = 0. 
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FIG. 11: The top row shows the real and imaginary parts (left and right, respectively) of the "doublet" of QNMs emerging 
from the algebraically special frequency as functions of a. The doublets only appear when rn > 0. We also ovcrplot the real 
and imaginary parts of the mode with I = 2, m = (showing the usual oscillatory behaviour). The bottom row shows, for 
completeness, the real and imaginary parts (left and right, respectively) of modes with negative m and branching from the 
algebraically special frequency. 



